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Photo-excited state in correlated electron system with spin- state degree of freedom is studied. 
We start from the two-orbital extended Hubbard model where energy difference between the two 
orbit als is introduced. Photo-excited met ast able state is examined based on the effective model 
Hamiltonian derived by the two-orbital Hubbard model. Spin-state change is induced by photo- 
irradiation in the low-spin band insulator near the phase boundary. High-spin state is stabilized 
by creating a ferromagnetic bound state with photo-doped hole carriers. An optical absorption 
occurs between the bonding and antibonding orbit als inside of the bound state. Time-evolution for 
photo-excited states is simulated in the time-dependent mean-field scheme. Pair-annihilations of 
the photo-doped electron and hole generate the high-spin state in a low-spin band insulator. We 
propose that this process is directly observed by the time-resolved photoemission experiments. 

PACS numbers: 78.20.Ls, 71.10.-w, 78.20.Bh, 78.47.J- 
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I. INTRODUCTION 

Optical properties and photo-induced phenomena in 
solids are one of the attractive themes in recent solid state 
physics. In particular, correlated electron system is one 
of the main targets for photo-induced exotic phenomena. 
Because of strong electron-electron interaction and multi- 
degrees of freedom, e.g. spin, charge, orbital and so on, 
a number of electronic and structural phases are realized 
under a subtle balance of interactions. 1 By irradiation 
of intensive laser pulse into one of the phases, a system 
is transferred into different phase transiently or perma- 
nently. This is termed photo-induced phase transition 
(PIPT) phenomena. 2 Nowadays, a number of experimen- 
tal and theoretical studies have been done in PIPT phe- 
nomena in transition- metal oxides^ - — low-dimensional 
organic salts j^ - — and others. 

Among the multi-degrees of freedom, spin-state degree 
of freedom has attracted much attention from view point 
of optical manipulation of magnetism. In a certain mag- 
netic ion, different magnitude of the spin angular momen- 
tum is realized by changing external fields, such as tem- 
perature, pressure and photons. This is termed the spin- 
state transition and is caused by a competition between 
the crystalline field splitting and the Hund's coupling. 
A well known example of the photo-induced spin state 
change is seen in the so-called spin cross-over complexes, 
such as Prussian blue analogue complex ji 2 - — Here, pho- 
tons bring about a charge transfer from the neighboring 
Fe ions to Co ions, associated with the spin-state change 
in Co ions from the low-spin (LS) state to the high-spin 
(HS) one. A main mechanism of the cooperative spin- 
state transition in a series of materials is supposed to be 
the elastic interaction^"— a local volume change of a 
metal-ligand cluster propagates over a crystal lattice. 

Another material where photo-induced spin- 
state change is realized is the perovskite cobaltites 
Ri- x A x CoO^ (R: a rear-earth ion, A: an alkaline-earth 



ion) and their families *22r— In an undoped compound 
LaCo03 j 23 i 24 a formal valence of a Co ion is 3+ with a 
d 6 electron configuration. There are possible three spin 
states: the LS state with the (t2 g ) 6 (e g )° configuration, 
the intermediate-spin (IS) state with (t2 g ) b (e g ) 1 , and the 
HS state with (t2 g ) 4 (e g ) 2 . It is supposed from the electric 
resistivity and the magnetic susceptibility measurements 
that the LS band insulator in low temperatures are 
changed into the HS or IS metallic state with increasing 
the temperature (T)*2§r— By substitution of R by A, 
corresponding to hole doping into the non-magnetic 
insulating ground state, a system shows ferromagnetic 
metallic behavior £t£l One key point to understand 
the electronic and magnetic properties in cobaltites 
are strong correlation between electron conduction and 
magnetism, i.e. charge and spin degrees of freedom of 
electrons. 33 

Optical irradiation and manipulation in perovskite 
cobaltites and related materials have been examined by 
the ultrafast optical pump-probe measurements j 34 i 35 Re- 
cently, detailed experiments and analyses have been done 
in so-called A-site ordered perovskite-type i?BaCo206.<5 
crystals by Okimoto and co-workers. 34 After pump pulse 
is introduced into the LS insulator, a metallic state, 
which is different from the high-temperature metallic 
state, is observed in the optical conductivity spectra. 
This photo-induced state strongly depends on the R 
species, which is supposed to control a ratio of electron 
correlation and band width. These experiments suggest 
that strong correlation between electronic and magnetic 
states remains even in the photo excited state, and tell us 
that the photo-irradiation phenomena in the cobalt ox- 
ides should be reexamined from different viewpoint from 
the photo-induced spin-state change in spin-cross over 
complexes. 

In this paper, photo-induced spin-state change in cor- 
related electron systems is studied theoretically. From 
the two-orbital Hubbard model, the effective Hamilto- 
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nian for the photo-excited state is derived. The photo- 
excited metastable state is obtained through analyses of 
the effective Hamiltonian by using the exact diagonal- 
ization method. By irradiation of photons into the LS 
band insulator near the phase boundary, the HS state is 
induced. It is found that the HS state is stabilized by 
forming a bound state with a photo-doped hole. This 
bound state brings about a characteristic peak structure 
in the optical spectra in the photo-excited state. A time 
evolution after photo- irradiation is examined in the time- 
dependent mean-field scheme. A creation of the HS state 
is caused by a pair annihilation of photo-doped electron 
and hole. This mechanism is able to be confirmed by the 
time-resolved photoemission spectroscopy experiments. 

In Sect. II, the model Hamiltonian and the effec- 
tive model for the photo-excited states are introduced. 
In Sect. Ill, numerical results of the electronic states 
before and after photo-irradiation are presented. In 
Sect. IV, the time-dependence of the photo-excited states 
are shown. Section V is devoted to discussion and con- 
cluding remarks. A brief report for the previous studies 
in the photo-induced metastable state was published in 
Ref . M 



II. MODEL 

A. two-orbital Hubbard model 

We start from the two-orbital Hubbard model as a 
minimal model to examine the photo-induced spin-state 
change. Two orbit als, termed A and B corresponding to 
the e g and t<i g orbitals in a Co ion, respectively, are intro- 
duced in each site in a lattice. The crystalline field split- 
ting between A and B is represented by A = ea~ > 
where sa and sb are the level energies of the A and B 
orbitals, respectively. The model Hamiltonian is given as 

U = Uu+U u (1) 

where we define the on-site term, 

i~f iaa' 
+ J C iAa C iBa' C iBa + I ^ C l^ C ^l Ci ^^ Ci ^ ' ^ 



and the inter-site term, 



(3) 



Here, c\ la is the electron creation operator at site i 
with orbital 7(= A, B) and spin <r(=t, \). We define 
the number operator rii 1CT = cj Ci 7<T and a subscript 
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FIG. 1: (color online) Local electronic configurations. 



7 = (A, B) for 7 = (B, A). The intra-orbital Coulomb 
interaction U ', the inter-orbital Coulomb interaction U' ', 
the pair-hopping /, and the Hund's coupling J are in- 
troduced. The electron transfer integrals between the 
nearest-neighboring (NN) sites are set to be diagonal 
with respect to the orbital. We assume a relation ts <^a 
by considering the transfer- integrals in perovskite oxides, 
and we take ^ = 1 as a unit of energy. 

Let us consider the local electronic structure, in which 
two electrons occupy each site, and the electron transfers 
are set to be zero. The following LS state with S = 
and the HS state with S = 1 are the possible ground 
states (see Figs. [TJa) and (b)). The eigen function and 
the eigen energy for the LS state are given as 



\*I>L) = (fAC ] M C ] Ai + /£<4 t <4 J |0) 



(4) 



and El = U + A — \/A 2 + J 2 , respectively, with the 
coefficients 



Sb 




-1/2 



and 



The wave functions for the HS state are given by 



I^h+i) = cUcU |0) , 



V2 

i)=cV f B jO>, 



(5) 

(6) 

(7) 
(8) 
(9) 



for S z =+1, and -1, respectively, and the energy in the 
HS triplet state is Eh = U' + A — J. In the numerical 
simulations, for simplicity, we assume the relations U — 
U' = 2 J, U = 4J and I = J. In these assumptions, the 
LS and HS states are degenerated at A = J/ \^3. 



B. Effective Hamiltonian 

One of the main purposes in this paper is to examine a 
stable steady photo-excited state. This is defined as the 
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FIG. 2: (color online) Adiabatic energy surfaces before and 
after photo-irradiation. Horizontal axis represents a number 
of the HS state. 



lowest-energy state inside of the energy surface, where 
density of the photo-excited electron-hole pairs is fixed. 
A schematic picture is shown in Fig. [2j Bold curves repre- 
sent the adiabatic energy surfaces before and after photo- 
irradiation as functions of a number of the HS sites. Pho- 
tons excite a system from the lowest-energy surface to 
the higher-energy surface. Through several kinds of re- 
laxation processes, the system settles down in the lowest 
energy state in the higher energy surface. Instead of time- 
dependent simulations for the photo-excited dynamics, 
we examine the lowest-energy state inside of the energy 
surface where a number of the electron-hole pair is one in 
a TV-site system. This state is termed the photo-induced 
metastable state, from now on. We derive the two ef- 
fective Hamiltonians, where numbers of the electron- hole 
pairs are zero and one in a TV-site system. 



Hamiltonian is given as 

H = Hu + Jhh ( S * ■ S i ~ X ) P ? P T 

m 

m 
m 

- J ++ VilJ (Si ■ Sj - 1) + (Si ■ S, - 1) 7+7+] 

{ij) 

+ J +" E i 1 ! • S i + !) T i + J r ( S i ■ Sj + 1) It] 



(10) 



Here, Si is the spin operator defined by Si = 
(1/2) J2 7 aa' c ly ff ^<r'V with the Pauli matrices a, and 
Pi and Pf 1 are the projection operators for the LS and 
HS state defined by 



and 



^ L = |fc)(fc| 



\*l>Hli) (i>Hli\ , 

l=(+,0,-) 



(11) 



(12) 



respectively. The operators 7+ and I i change the spin 
state as 



Li\ 



and 



(13) 



(14) 



The prefactors in each term in Eq. ([TP]) are the exchange 
constants defined by 



J ™ = -WTJ ' 



Jll = 



2U r + 2Aj - U - J' 



(15) 



(16) 



U' + Aj- A /7 / + A J + A J ' 
The effective Hamiltonians are derived by the pertur- (17) 
bational processes from the two-orbital Hubbard model 
in Eq. (pQ). The inter-site transfer term, T-L t , is treated as 

the perturbation term. As for the effective Hamiltonian j ++ _ 2tAtBfB Ia 



before the photo-irradiation, the HS and LS states, de- \U + J 2U' — U — J + 2 A j J 
fined in Eq. ^ and Eqs. ©-(|9]), respectively, are adopted (18) 
as the basis states. Other local-states, where two elec- 
trons occupy each site, have higher energies of the order . 2 2 . 

of A, J, and U than the LS and HS sates. By consider- = 2tAtB ( I ) > (19) 

ing all of the second-order perturbational processes, the \U f + Aj — A U' + Aj + Ay 



where we define Aj = ^ A 2 + J 2 . 

The effective Hamiltonian after photo- irradiation is de- 
rived in the same way. As the unperturbed states, in ad- 
dition to the LS and HS states, we introduce the states 
where numbers of electrons in a site are one or three (see 
Figs. QJc) and (d)). These local states are termed the 
hole state and the electron state, respectively. The wave 
functions are given as 



10} , 



and 



(20) 



(21) 



respectively. The eigen energies are E e = A-\-U + 2U f — J 
for the electron state and Eh = for the hole state. We 
assume that a number of both the electron state and 
the hole state is one in a TV-site cluster. The calculated 
effective Hamiltonian is classified by the electronic states 
in the NN sites as 



Hi — ^0 + T~ieh + H e + T~ih- 



(22) 



The first term Ho corresponds to Ho in Eq. ([TP]) , where 
both the electron and hole states are not concerned in 
the interactions. The second term is for the interactions 
between the electron state and the hole state. The third 
and fourth terms describe the interactions between the 
electron state and LS or HS, and the interactions between 
the hole state and LS or HS, respectively. Explicit forms 
for the Hamiltonian are given in Appendix. 

The ground state before photoirradiation and the 
photo-induced metastable state are obtained in the ef- 
fective Hamiltonians in Eqs. ([T0|) and (|22]h respectively, 
which are analyzed by the exact-diagonalization method 
based on the Lanczos algorithm. Time evolutions in the 
photo-induced dynamics are calculated in the two-orbital 
Hubbard model in Eq. (pQ). 



III. ELECTRONIC STATES BEFORE AND 
AFTER PHOTOIRRADIATION 



CO 



1 

0.8 

0.6 

0.4 

0.2 


0.6 

0.4 

0.2 V 



0.5 



(a) 



LS 


LS 


LS 


LS 


LS 


LS 


LS 


LS 


LS 



HS 


HS 


HS 


HS 


HS 


HS 


HS 


HS 


HS 



LS 


HS 


LS 


HS 


LS 


HS 


LS 


HS 


LS 



(b) 



Q 

(0,0) 
(tt,0) 

(7T,7r) 



(c) 




3.2 



3.4 3.6 

J/t A 



3.8 



FIG. 3: (color online) (a) The number density of the HS 
state riHS, (b) the spin correlation function S(q), and (c) the 
spin-state correlation function I(q) in the ground state. The 
parameter values are chosen to be U — 4 J, U' — 2 J, A = 
lOt a, and tB — 0.05£a- A two dimensional cluster of the 
N = 8 sites with the periodic boundary condition is adopted. 



A. Ground State 

Electronic structure in the ground state is examined 
by analyzing the effective Hamiltonian Ho in a finite 
size cluster system. Several physical quantities are plot- 
ted in Fig. [3] as a function of the Hund's coupling J at 
A/tA = 10. We introduce the number density of the HS 
states which are estimated from the electron number in 
the orbital A defined by 

the spin correlation function, 

S ^ = ^E e-^~^ {S, -Si), (24) 



and the spin-state correlation function defined by 

J (9) = ^E e "' 9 ' ( "" , ' 3) ( J ^>- ( 25 ) 

Here, we define the spin-state operator as a projection 
operator by 

I! = \ E (\TpHmi)(lpHmi\-\lpLi)(lpLi\), (26) 
ra=±l,0 

which takes 1/2 and —1/2 for the HS and LS states, 
respectively. With increasing J, three different phases 
appear in Fig. [3j In a region for small J, both uhs an d 
S(q) are zero, and 7(0, 0) is almost one. On the other 
hand, in a region of large J, nns and 7(0, 0) are one, and 
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FIG. 4: (color online) Phase diagram in the plane of the crys- 
talline field splitting A and the Hund's coupling J. Broken 
and bold lines represent the phase boundaries in the ground 
state and in the photo-induced metastable state, respectively. 
Abbreviations, HS, LS, and H/L represent the HS phase, the 
LS phase, and the HS-LS mixed phase, respectively. A verti- 
cal dotted line represents the parameter region where the data 
in Figs. [3] and [5] are calculated. Parameter values are chosen 
to be U = 4 J, U' = 2 J, and t B = 0.05£a- A two dimen- 
sional cluster of the N = 8 sites with the periodic boundary 
condition is adopted. 




5(7r, 7r) is the largest. Two phases are identified as the 
LS band insulator and the HS antiferromagnetic Mott in- 
sulator. Between the two, there is an intermediate phase 
where uhs — 0.5, and 7(7r,7r) is one. These data imply 
that the HS and LS states are aligned alternately. This 
phase is termed the spin-state ordered phased This al- 
ternate ordering of the HS and LS states is caused by the 
fourth term in the right hand side in Eq. ([TO]) ; Jlh given 
in Eq. (fT7|) represents the attractive interaction between 
the LS and HS states. 

Numerical data in several J and A are summarized 
in a phase diagram shown in Fig. 2J where the phase 
boundaries in the ground state and those in the photo- 
induced metastable state are plotted in one figure. The 
detailed results in the metastable state will be presented 
in Sec. IIIIBI Here, we identify the LS (HS) phase as a 
state, where electron numbers of the A orbital is smaller 
(larger) than 0.3. We confirm that the size dependence 
of the phase boundaries is of the order of 0.01^. The 
LS and HS phases appear in regions of large A and large 
J, respectively. The spin-state ordered phase appears 
between the two phases. 



FIG. 5: (color online) (a) The number density of the HS 
states riHS, (b) the spin correlation function S(q), and (c) 
the spin-state correlation function I(q) in the photo-excited 
metastable state. In comparison, results of nus in the ground 
state are also plotted in (a). Parameter values are chosen to 
be U = 4 J, U' = 2J, A = lOt a, and t B 0.0bt A - A two 
dimensional cluster of the N = 8 with the periodic boundary 
condition sites is adopted. 



B. Photo-induced Metastable State 

Several physical quantities in the photo-induced 
metastable state are presented in Fig. [5] as a function of 
the Hund's coupling J at A/tA — 10. The number den- 
sity of HS states is estimated from a number of electrons 
in the orbital A defined by 

nm = (tha) - 1 J , (27) 

where the electron and hole states are subtracted in a 
denominator. It is shown in the Fig. [5fa) that nHS 
in the photo-induced metastable state is finite between 
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3.25 < J/tA < 3.30 where uhs m the ground state is 
zero. Different value between the two states implies that 
one HS state is generated in the TV-site cluster. This 
phase in the photo-induced metastable state is distinct 
from the spin-state ordered phase observed in the ground 
state; the spin-state correlation functions at any q are not 
remarkable, and a weak spin correlation at q = (0, 0) is 
observed. Detail properties of this phase are introduced 
latter. 

The phase diagram in the photo-excited metastable 
state is presented in Fig. HJ together with that in the 
ground state. The phase boundary between the LS and 
LS-HS mixed phases shifts to a region of the LS phase. 
There is a parameter region where the LS phase in the 
ground state is changed into the LS-HS mixed phase in 
the photo-excited metastable state. That is to say, the 
photo-irradiation induces the HS state in the LS phase 
at vicinity of the phase boundary. We note that the spin- 
state change also occurs from the HS phase to the mixed 
phase. 

Now we examine the electronic structure in the photo- 
induced HS state in more detail. We introduce the 
electronic-state distribution function defined by 



9M(n) 



jenNN i 



I pM p h\ 



(28) 



where J2j enNN implies a summation of j connecting the 
n-th NN sites of z, and z n is a number of the n-th NN 
sites. The operator (M = L, H, e) is the projection 
operator for the M-state at site i. The operators for 
the LS and HS states are define in Eqs. ([TT]) and ([T2]) . 
respectively, and those for the electron and hole states 
are defined as 



and 



P t k = Yl \^ai) Who 



(29) 



(30) 



(a) 



J, X, q X, .-, 

x xT^ x 
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FIG. 6: (color online) (a) A schematic definition of the dis- 
tribution function gM(n). (b) Distribution function of the 
electronic states as functions of a distance from the hole state. 
Bold, broken and dotted lines represent the distribution func- 
tions for the HS, LS and electron states, respectively. Two 
dimensional N = 10 site cluster with the periodic boundary 
condition is adopted. The inset shows distribution functions 
for the HS state in one dimensional N =6, 8, and 10 site clus- 
ters with the periodic boundary condition. Parameter values 
are chosen to be J = 3.3£a, U = 4 J, U' — 2 J, A = 10£a, and 

t B = omt A . 



respectively. This function, ^m(^), describes distribu- 
tion of the local electronic states at the n-th NN sites 
from the photo-induced hole state, as shown in Fig.[6ja). 
Numerical results of the distribution functions in a two- 
dimensional cluster are shown in Fig. [6fb). Parameter 
values are chosen to be J = 3Ma and A = 3Ma in which 
the HS state is induced by photo-irradiation. A charac- 
teristic feature is shown in the HS distribution function; 
9H{n) is nearly 0.25 at n = 1 and zero at n > 2. This 
implies a local bound state between the HS state and a 
photo-doped hole state. The size dependence of ^(n) is 
checked in the one-dimensional clusters, and results are 
shown in the inset of Fig.EJb). The HS distribution is lo- 
cated at the NN sites of the hole state, and almost no size 
dependence is seen in the results. Different numerical val- 
ues of ^(n) in one and two dimensional clusters, i.e. 0.25 
and 0.5, are attributed to difference of z n . Spin structure 



in this bound state is monitored by the correlation func- 
tion defined by Yli^j {&j ' ^i^j 1 ^i) which represents 
the spin correlation between the hole and HS states. Cal- 
culated value is about 0.5 which implies a ferromagnetic 
spin correlation. Figure[6]also shows that g e (n) monoton- 
ically increases with n. This is due to the kinetic-energy 
gain of the photo-excited electron. A schematic electronic 
structure in the photo-induced metastable state is pre- 
sented in Fig. [71 

Here we discuss a mechanism of the ferromagnetic HS- 
hole bound state. In the ground state, the energy dif- 
ference per site between the LS state and the HS state 
is given by AE HS - ls = E HS - E LS = (U' + A - J) - 
(U + A — VA 2 + I 2 ) in the local limit. Let us consider a 
situation that photo-excited electrons and holes are intro- 
duced in the LS phase at vicinity of the phase boundary, 
and these photo-carriers move in the system. Magnitude 
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FIG. 7: (color online) (left) A schematic picture of the HS- 
hole bound state, (right) Energy levels inside of the HS-hole 
bound state. 




agram. In FiglHJ the phase diagrams in the ground 
state and the photo-induced metastable state are plot- 
ted as functions of a ratio of the band widths for the A 
and B bands, i.e. t B /tA- In a region of tsj^A < 0.1, 
there is a phase space where the LS phase in the ground 
state is changed into the HS-LS mixed phase in the 
photo-induced metastable state. With increasing tsl^A-, 
this phase space is shrunken and disappears. This is 
explained from above consideration where stability of 
the photo-induced HS state is controlled by a factor 
AEbs — {~^a) — {—zf^ts)- This tendency for stability 
of the HS state is similar to the previous results in the 
chemical doping^ 



C. Optical Spectra 




0.1 
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FIG. 8: (color online) Phase diagrams in the ground state 
and the photo-excited metastable state in the plane of tB/tA 
and J It a- A shaded area implies a parameter region where 
the HS state is induced by photo-irradiation. Two dimen- 
sional N = S site clusters with the periodic boundary condi- 
tion is adopted. Parameter values are chosen to be U = 4 J, 
U' 2J, and A 10t A . 



of the exchange process between the hole state and the 
LS state is given by f B ^B (see Eq. (|A2|) ). Thus, when the 
electron and hole states move around the LS background 
without generation of the HS state, the kinetic energy of 
the hole state is —zf B tB where z is a number of the NN 
sites. On the other hand, when the ferromagnetic HS- 
hole bound state is generated, one electron occupies the 
bonding orbital in the bound state and this energy gain 
is — Ia as shown in Fig. [7J Since the electron states move 
in the same ways in both the two cases, energy differ- 
ence are AEbs = (— ^a) — {—zf B t B ). When this energy 
gain due to the bound state, \AEbs\, overcomes the en- 
ergy cost for the HS generation, AEhs-ls, the HS-hole 
bound state is realized. 

Above consideration for the energy balance is con- 
firmed in the band width dependence of the phase di- 



In this subsection, we show the optical spectra in the 
photo-induced metastable state. The optical absorption 
spectra are defined by 



a 



-—im^oir- 

Ntt uo ■ 



H e ff + E 



IT] 



f |Vo> , 

(31) 



where is taken to be Ho in Eq. (fTD|) and Hi in 

Eq. ([22]) for the ground state and the photo-induced 
metastable state, respectively, |^o) and Eq are the cor- 
responding lowest energy state and energy, respectively, 
and a(f3) represents a Cartesian coordinate. We intro- 
duce the current operator 



7 — ^ ^ ^ 7 (c\^ a Ci+aya H.C.^ 



(32) 



which is defined in the restricted Hilbert space in each ef- 
fective Hamiltonian. A damping constant is introduced 
as f]. The optical spectra are calculated by the exact 
diagonalization method based on the recursion proce- 
dure. Two dimensional finite-size clusters with the peri- 
odic boundary condition is adopted. 

The absorption spectra in the photo-induced 
metastable state, where the HS state is induced by 
photo-irradiation, are shown in Fig. [9] The system size 
is taken to be N = 8 and 10. Characteristic two peaks 
appear in the spectra at uj = 2. It a and uo = 2Ma in the 
case of N = 8. These are termed the peaks B and D, 
and their energies are denoted ujb and ujd, respectively. 
In order to assign these peaks, we calculate the bond 
correlation function in the excited states given as 



E 

(r/>7<7 



( pm pi J pi prr 



H.c.) |VK)>. 

(33) 



where \^(uj n )} is the eigen function of the Hamiltonian 
Hi corresponding to the final state of the n-th optical ab- 
sorption peak, and uo n is its eigen energy. Eigen functions 
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FIG. 9: (color online) Optical absorption spectra in the 
photo-induced metastable state, where the HS is induced. 
Bond correlation function B^ l,rn \oj) is also plotted. Red bro- 
ken lines and blue dotted lines are for (Z,m) = (h,HS) and 
(/, m) — (e, LS), respectively. Cluster size is N — 8 in (a), and 
N — 10 in (b). Parameter values are chosen to be J — 3.3£a, 
U = 4J, U' = 2J, A = lO.OtA, t B = 0.05tA, and r] = 0.2t A . 



and eigen energies are obtained by using the conjugate 
gradient method. This function measures the bond cor- 
relation between the I and m local electronic states in the 
photo-excited state. 

Numerical results of this correlation function together 
with the optical absorption spectra are presented in 
Fig. [91(a), where we set (Z,m) = (h,HS) and (e,LS). In 
the ground state, i.e. uj = 0, B^ h,HS \uj = 0) ~ —1 and 
B^ e ' LS \uj = 0) ~ —3. These values are consistent with 
the picture presented in Fig. [7J where a photo-doped hole 
forms a bound state with HS, and a photo-doped electron 
is located in a bottom of the A-orbital band. In the ex- 
cited state corresponding to the peak B, B^ ls \ujb) — 
B^ LS \0) and B^ hs \lj b ) ~ 1 > B^ HS \0). This 
value of B^ h ' HS ^ (ujb) is interpreted that an electron occu- 
pies the antibonding orbital in the HS-hole bound state, 
and the peak B is assigned as an excitation between the 
bonding and antibonding orbitals inside of the bound 
state. As for the peak D, B^ ls \uj d ) > B^ LS \0) and 




FIG. 10: (color online) Optical absorption spectra in the 
photo-induced metastable state, where HS state is not in- 
duced. Broken lines represent the spectra obtained by the 
HC fermion model. Cluster sizes are N = 8 and 10. Inset 
shows size dependences of the peak energies in the optical 
absorption spectra calculated in the HC fermion model. Pa- 
rameter values are chosen to be J = 3.1t A , U = 4J, U' = 2J, 
A = 10.0£a, ts — 0.05£a, and 77 = 0.2£a in the original model, 
and t B = 0mtA,t ex = 10" 10 tA,ry = 0.2t A in the HC fermion 
model. 

B( h > HS )(uj D ) ~ B^ HS \0) which imply that a change 
in the photo-doped electron motion is concerned in this 
peak. 

Results in different size cluster of N = 10 are shown 
in Fig. EJb). Numerical values of B^ HS "> and B^ LS "> 
are almost same with the values in N = 8. Energy of 
the peak B is almost unchanged, but that of the peak 
D decreases with increasing N . These size dependences 
are consistent with the assignments that the peak B is 
attributed to the local excitation, and the peak D is re- 
lated to the kinetic motion of the photo-excited electron. 
We further examine the size dependence of the peak po- 
sitions in the one-dimensional clusters, and observe that 
the energy of the peak D decreases with the system size^£ 
This peak is interpreted as a Drude-like component in the 
thermodynamic limit. 

The optical absorption spectra in the photo-induced 
metastable state, where HS is not induced, are presented 
in Fig. [10] (see bold lines). Two peak structure is observed 
in N = 10 and are almost overlapped around uo = 3. Ha 
in TV = 8. Results are well reproduced by the following 
hard-core two-fermion model defined by 

Uhc = ~t A ^K f %- + H - c -) - tB ^2(b\bj + H.c.) 

-t e£C ^(aJ&]Mi+H.c), (34) 

(ij) 

where ai and fr^are the spin- less fermion operators at site 
i and describe annihilations of the electron state and the 
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hole state, respectively. We take a condition of a\b\ = 0. 
The first and second terms represent kinetic motions of 
the electron and hole states in the LS phase, respectively, 
and the third term represents an exchange of the electron 
and hole states. This model is derived in the limiting 
case of A >> / as follows. From Eqs. (|4])-(j6]), we have 
= c Bt c B| |0) wm ch is set to be a vacuum, |0), in this 
model. The electron and hole states are defined from this 
vacuum as \e) = a\ |5) and \h) = b \ respectively. The 
matrix elements for the exchange of the electron (hole) 
and LS states, corresponding to the first (second) term 
in Eq. (|3l]), are give n by t A j\ ~ t A (t B fg ~ t B ) from 
Eqs. (|A2|) and (|A12[) . The exchange of the electron and 
hole states, corresponding to the last term in Eq. (|34|h 
are given in the matrix elements in Eqs. (jA6|) and (|A9|) . 
We confirm numerically that this contribution to the op- 
tical spectra is much smaller than other terms, and set 
to be a small constant t ex = 10 _10 £a in the numerical 
calculation. In this effective model, the current operator 
along an a direction is given by 

Jhc = il A ^2(a\a i+a - H.c.) - it B ^2(b\b i+a - H.c), 

i i 

(35) 

The optical absorption spectra obtained in the HC 
model are shown by broken lines in Fig. [TUJ One pair 
of the a and b fermions is introduced in the N-site clus- 
ters. Spectra in the effective Hamiltonian are well repro- 
duced by the HC model. Size dependences of the peak 
positions are examined in detail in this model. Two di- 
mensional clusters with N = 8, 10, 4 x 4, 6 x 6, 8 x 8, 
10 x 10, 12 x 12 with the periodic boundary condition 
are adopted. Peak energies for the lowest three peaks 
are plotted in the inset of Fig. [10] as functions of N -1 . 
Energies tend to be zero in the thermodynamic limit. We 
interpret that these peaks originate from the metallic be- 
haviors of photo-doped electron and hole. 



IV. TIME-DEPENDENCE OF 
PHOTO-EXCITED STATE 



In this section, we show the real-time evolution 
of photo-exited state calculated in the mean-field 
scheme ^ 38 ! 39 and reveal a mechanism of the photo- 
induced HS state. 



A. Formulation 

Time evolution of the photo-excited state is analyzed. 
A mean- field type decoupling is applied into the Coulomb 
and exchange interaction terms in the two-orbital Hub- 



bard model in Eq. (pQ) as follows, 

i'ya L cr' 

( n nt) ( n ill) + U ' ^ ( n iAa) (n iB cT<) 

+ H.c.) 

icr (ij)l<7 

+ A^n m , (36) 



where (• • • ) implies the average calculated by the time- 
dependent mean-field wave function, and a subscript a 
is defined by a = (t,|) for a = (|,t). We note that the 
pair-hopping interaction in the Hamiltonian and the Fock 
terms are not taken into account. This is essential to re- 
produce the electronic states in the case of t A = t B = 0. 
The initial electronic wave function before the photo- 
excitation is obtained by solving the self-consistent equa- 
tions. The photo- irradiation is simulated by excitations 
of electrons from the highest occupied levels to the low- 
est unoccupied ones at time r = with conserving the 
^-component of the total spin-angular momentum and 
the total momentum. The time evolution of the wave 
function is calculated in the time-dependent mean-field 
scheme. The time-dependent Shrodinger equation for the 
u-th. level, \(j) u (r)), is given as 



d T "H M F{T') 



1^(0)), (37) 



where HmfW is the time-dependent Hamiltonian given 
in Eq. ([36]) . and P is the time-ordering operator. The 
wave function at time r + dr, where dr is short time 
distance, is calculated from the wave function at time r 
by expanding the exponential factor as 

\Mr + dr)) = (Mr)\Mr)) e~^ dT |^(r)> , 

/I 

(38) 

where \(p^(r)) is the eigen state of Hmf(t) with the 
eigen energy s^{t). In the numerical calculation, we take 
drt A = 10 -3 ~ 10 -4 , and check that the total energy is 
conserved within the order of 10 -2 percent. 

Phase diagram in the ground state is presented in 
Fig. [TTJ Phase boundaries are determined by the HS 
density. We also plot the results obtained by the exact 
diagonalization method applied to the effective Hamil- 
tonian shown in Fig. |4j Two results are qualitatively 
similar with each other, although the LS phase in the 
present calculation shifts to a low J region. 

B. Numerical Results 

Time evolution of the photo-excited electronic state is 
examined in a two dimensional N = 10 x 10 site cluster 
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FIG. 11: (color online) Phase diagram in the ground state 
obtained by the mean- field approximation. Two dimensional 
N = 10 x 10 site cluster with the periodic boundary condition 
is adopted. As a comparison, phase boundaries obtained by 
the exact diagonalization method on the effective Hamiltonian 
(N — 8) are plotted by broken lines. Abbreviations, HS, LS, 
and H/L, represent the HS phase, the LS phase, and the HS- 
LS mixed phase, respectively. Parameter values are chosen to 
be U = 47, U' = 2J, and t B = 0.05U- 



with the periodic boundary condition. A number of pho- 
ton in the cluster is chosen to be 7V p h = 10, which are 
introduced into the LS phase at a vicinity of the bound- 
ary (J = 3.Ha, A = 10t^). We monitor numbers of the 
HS state, the photo-doped electron state, and the photo- 
doped hole state, by the following physical quantities, 



with 



Ng\ iV e = and N h = £ iCT ^ 



Nf a = {riiAa) (1 - {niAa)) {niBa) {riiBa} , 



(39) 



Nl = (1 - (n iAa ))(l - {n iA - a )) (n iBa ) (1 



(niBa)), 

(40) 



and 



HS 



(niAa} (1 - (niAa)) (^Ba) (1 ~ (niBa)), (41) 



respectively. It is noticed that these are defined in prod- 
ucts of the mean-field number density in each orbital 
and spin, instead of the projection operators such as Pf 
(Eq. (|29|)), if (Eq. d30j) ) and P % H (Eq. (H2J)), which can- 
not be calculated directly in the mean- field scheme. 

Time dependence of these numbers are plotted in 
Fig. [T2l a). Except for the early time below rtA = 10, 
where all three are almost constant, N HS increases, and 
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FIG. 12: (color online) (a) Time evolutions of numbers of 
the HS state, the electron state and the hole state, (b) A 
schematic picture of the electron-hole pair annihilation pro- 
cesses. Two dimensional N = 10 X 10 site cluster with the 
periodic boundary condition is adopted. Parameter values 
are chosen to be U = 4 J, U' — 2J, J — 3. It a, A = 10£a, 
t B = 0.05£a, and N ph = 10. 



N e and N h decrease monotonically. That is, changes in 
the three numbers are correlated with each other. This 
result is interpreted that the HS states are created by an- 
nihilation of the photo-induced electron and hole states. 
Let us consider a situation where the electron and hole 
states adjoin, as shown in Fig. fl2T b). When an electron 
in the A orbital transfers to the hole state, a LS-HS pair 
is generated. This pair is also generated by the electron 
transfer in the B orbital. This is termed the electron-hole 
pair annihilation process, from now on. 

Snapshots for the local electronic states are shown in 
Fig. [T3j At r = lOt^ 1 , three numbers are almost ho- 
mogeneous. At r = lOOt^ 1 , distributions for the electron 
states start to be inhomogeneous and a vertical shape do- 
main appears in N^. At r = SOO^ 1 , the electron states 
in this vertical-shape domain begin to be localized. It is 
shown that, at the sites where is large, and N HS 
are also large. These data support the mechanism of the 
HS creation due to the electron-hole annihilation pro- 

at 



cess. After the HS states are generated, and iVJ 



t 

the same sites still remain to be large. This observation 
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rt A = 10 rt A - 100 rt A - 500 



rt A - 900 




FIG. 13: (color online) Snapshots of the electron state, the 



hole state and the HS state 

-1 irnnj.-! j 



Time is chosen to be lOt , 



lOOt^ 1 , 500£^ and 900£^\ Two dimensional N = 10 x 10 
site cluster with the periodic boundary condition is adopted. 
Parameter values are chosen to be U = 4J, U' — 2J, J — 
3.1t A , A = 10£a, t B = 0mt A , and iV ph = 10. 



is not contradict to the electron-hole annihilation, but is 
due to the definitions of and [see Eqs. (|5Dj) and 
([39]) ]; these are represented in the products of the mean- 
field number density, instead of the projection operators. 
We have numerically confirmed that, as Nf and Nj 1 in- 
crease, N^ s increases, and then Nf and N- 1 decrease. 
Relation between the numerical results of the real-space 
snapshots and the HS-hole bound state introduced in the 
previous section will be discussed in Sect. El 

The electron-hole pair annihilation processes are also 
examined by the time dependent density of state (DOS). 
We define DOS as 



A(u) = A e {oj)+A h {uj), 
with the electron part 

A%co) = Y / ^-eAr)} (etc,), 

and the hole part 

A h {oo) = Y J ^-^{r)] (c u ct). 



(42) 



(43) 



(44) 
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FIG. 14: (color online) Density of states. The electron 
part and the hole part of DOS are represented by pink solid 
and blue solid lines, respectively. Time is taken to be lOt^ 1 , 
lOOt^ 1 , 500t^\ and 900^ . Two dimensional N = 10 x 10 
site cluster with the periodic boundary condition is adopted. 
Parameter values are chosen to be [/ = 4J, U' — 2J, 
J 3.1£a, A 10£a, t B 0.05U, r\ 0.1t A , and N ph 10. 



The operator cj, is the creation operator obtained by di- 
agonalizing the Hamiltonian at time r, e v {r) is the corre- 
sponding mean-field energy, and (• • • ) implies the average 
in terms of the wave function of \ip(r)). In the numerical 
calculation, the delta functions in Eqs. (|43|) and ([44]) are 
replaced by the Lorentz function with a damping con- 
stant T] = O.Ha- 

Numerical results of the time- dependent DOS are 
shown in Fig. [14j At r = lOt^ 1 , an energy gap exists 
between the narrow B band and the wide A band. Tiny 



weights of the hole and electron parts of DOS are ob- 
served in the top of the B band and the bottom of the 
A band, respectively. At time r = 500^, the top of the 
B band and the bottom of the A band start to separate 
from the main bands. Finally, at r = OOOt^ 1 , the original 
gap is almost filled out by in-gap states. These data are 
consistent with the localization of the electron and hole 
states observed in the snapshots in Fig. [131 

Based on these results, we consider the energy balance 
in the electron- hole pair annihilation processes. On-site 
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mean-field energies of the LS, HS, electron, and hole 
states are given by E^p = U, E^§ = U' + A - J, 
^mf = U + 2U' + A - J, and = 0, respectively. 

When one electron-hole pair is changed into one LS 
state and one HS state, the on-site energy is changed 
as AE eh ^ LH = (E{& + E*§) - (E^ F + ^ F ) = -*7' 
which is negative, i.e. energy loss. This energy is com- 
pensated by the kinetic energy of the hole and electron 
states, which are not concerned in the pair-annihilation 
processes. This is confirmed in DOS at rt^ = 900 (see 
Fig. H3Td)): the hole part of DOS in the B-orbital band 
distributes not only to the top of the band, but also down 
to the middle of the band. This indicates increasing of 
the kinetic energy of holes with time. 

In the last part of this section, we examine, on 
the time-evolution of the photo-induced HS generation, 
roles of the relativist ic spin-orbit (SO) interaction which 
breaks the spin angular-momentum conservation. Here 
we mimic the SO interaction in the 3d orbit als as follows, 

ri S =i£j2 (°\a^H + 4ai C ibt " 4^*4 " 4| C ^t) > 

(45) 

with the SO interaction constant £. It is demonstrated 
that, when this interaction acts on the LS state, the HS 
state is created as follows, 

Kso4 bt 4 H |0> = -i£ (4 t 4 t - 4ai4 H ) i°) • ( 46 ) 

Numerical results of the time evolutions of N HS and 
N h are presented in Fig. [15j where the SO interaction 
constant is taken to be £ = and OMa- Before r = 
200£^ 1 , the SO interaction effects are not seen in N HS . 
However, beyond r = 200£^ 1 , N HS starts to decrease in 
the case of a finite £. The observed reduction of the HS 
state in the case of finite £ is due to the transition from 
the HS to LS states through the SO interaction, as shown 
schematically in Fig. fT5l b). This result indicate that 
roles of the SO interaction on the spin-state transition is 
destructive rather than constructive. 



V. DISCUSSION AND CONCLUSION 

In this section, we remark i) a connection between the 
calculated results in the photo-excited metastable state 
and the time dependent simulation in the photo-excited 
state, which are presented in Sects. Hill and [IVl respec- 
tively, and ii) implications of the present theoretical re- 
sults to the recent experiments. 

In Sect. II VI we show in the time-dependent simulation 
that a pair annihilation of photo-doped electron and hole 
generates a HS state. In this scheme, here we discuss a 
stability of the HS state and a role of the photo-doped 
hole. A correlation between the local HS state and the 
hole state around the local HS state is examined numer- 
ically. We introduce the number density of the local HS 




FIG. 15: (color online) (a) Numbers of the HS state, the 
electron state and the hole state in the model with the spin- 
orbit interaction. Two dimensional N = 10 x 10 site cluster 
with the periodic boundary condition is adopted. Parame- 
ter values are chosen to be U = 4 J, U' = 2 J, J = 3. It a, 
A = lOt a, ts = 0.05£a, and N p h = 10. Ten data sets with 
different initial values for the time evolutions are averaged, 
(b) A schematic picture of the transition between the HS state 
and the LS state due to the SO interaction. 



state at time r, N^ s (t), defined in Eq. (j4TJ), and the 
number density of the hole state around the site i at time 
t' denoted as p ia (r f ) = E^jU^'). We define A^(r') 
in Eq. (|^0|) . and introduce a symbol Ej which implies a 
summation for the NN sites of i. Data sets of ^ a (r) and 
iV^ 5 (r + Ar) with At = lOt^ 1 are obtained in 100 times 
simulations with different initial states. Numerical re- 
sults are shown in Fig.[16j A positive correlation between 
the two quantities is seen in this figure. In particular, in 
a region of high density of hole, {p% a {j) ^ 0-1)? a number 
of data for small N^ s (r + At) is a few. On the other 
hand, in a region of low density of hole, (picr(r) < 0.1), 
value of N^ s (t+At) distributes. These results are inter- 
preted that in the case that the hole density around the 
photo-induced HS state is low, a probability of the sur- 
vival HS state number is randomly distributed. On the 
other hand, in the case of high hole density around the 
HS state, the HS density increases with increasing the 
hole density. These relations between the photo-doped 
hole and the HS state are consistent with the results in 
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A hole 




r > t + At 



FIG. 16: (color online) Correlation between the number 
density of the local HS state at time r, N^ s (r + At), and 
the number density of the hole state around the site i at 
time r', pia(r). We take Ar = 10/£a- Two dimensional 
N — 10 x 10 site cluster with the periodic boundary condi- 
tion is adopted. Parameter values are chosen to be U — 4 J, 
U' — 2J, J — 3.1£a, A 10£a, t B 0.05£a, N ph 10, and 
Ar = lOt^ 1 . 

Sect. IHIi where the HS state is stabilized by forming the 
HS-hole bound state. 

Next we compare the present calculated results with 
the experimental data reported in Ref. [34|. As introduced 
previously, key points in the optical pump-probe exper- 
iments in i£BaCo206.5 are a) a photo-induced metal- 
lic state is different from the high-temperature metallic 
state, and b) this photo-induced state strongly depends 
on the R species. From the calculated results, we propose 
that the observed metallic state is attributed to the HS- 
hole bound state. Experimental spectral weight induced 
by the photon pumping is interpreted to be the dipole 
transition inside of the bound state. We have checked 
by the exact diagonalization method in a small size clus- 
ter that a clear bound state between thermal hole carri- 
ers and the HS states is not stabilized in finite tempera- 
tures^ It is well known that, in the perovskite crystal, 
the electron transfer intensity is systematically controlled 
by the R species through a changing of a Co-O-Co bond 
angle. The smaller the ionic radius of the R ion is, the 



smaller the e g band width is. With increasing the ionic 
radius from Tb to Sm, the photo-induced metallic state is 
remarkably seen in the experimental optical conductivity 
spectra. These data correspond to the calculated results 
in the phase diagram in Fig. |4j Increasing of the transfer 
integral of the A band, being equivalent to decreasing of 
A/tA, is indicated by the arrow in this phase diagram. 
A system is transferred from the phase, where the spin- 
state is not changed by photo-excitation, to the phase, 
where the HS state is induced by photo-irradiation. This 
consistency between the theory and the experiments is 
additional evidence of existence of the photo-induced HS- 
hole bound state. 

In conclusion, we study the photo-induced spin-state 
change in a correlated electron system. The photo- 
induced metastable state is examined in the effective 
Hamiltonian which is derived by the two-orbital Hubbard 
model. By photo-irradiation into the LS phase near the 
phase boundary with the mixed phase, the HS state is 
induced and is stabilized by forming a bound state with 
a photo-doped hole state. The optical transition inside 
of this bound state appears. Time dependent simulation 
for the photo-excited state is also performed on the two- 
orbital Hubbard model in the time-dependent mean-field 
scheme. A pair annihilation of the photo-doped electron 
and hole states generates the HS state. This process re- 
flects on the time-dependent DOS. The present results 
propose a new state of the photo-excited matter in cor- 
related electron system with multi-degrees of freedom. 
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Appendix A: Effective Hamiltonian for the 
Photo-excited Meta-stable State 

In this Appendix, explicit formulae of the effective 
Hamiltonian for the photo-induced metastable state are 
presented. Matrix elements in terms of the electronic 
states in NN sites are shown. These are classified by the 
electron number, n, and the z-component of the total 
spin-angular momentum, S z , as l-L^ n ' Sz \ The wave func- 
tions in the two sites are denoted as In the fol- 
lowing notation, each term in the Hamiltonian in Eq. ([22]) 
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is given as U eh = H (4,0) , U e = H (5,3/2) +^( 5 ' 1 / 2 ) 

and 74=^ (3 ' 3/2) +^ (3 ' 1/2) - 
(1) (n = 3, 5, = 3/2) 



^(3,3/2) 



*A 
*A 



The basis set is {\^^ H+1 ) \^ H+1 ^ H ^)}. 
(2) (n = 3, S z = 1/2) 
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where J eh ± = J ehl ± J e/l2 and 



JWi2 — (^a + 
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with 



where 




-1/2 



(A5) 
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Oieh 



'- A- A J' 



(A7) 



(A8) 



The basis set is {|^ et ,^ t ) l^t^et)}- 



and 



t 2 f 

l AJA 



2Aj 

JA l B 



A + 2C/' + J + Aj 

3 



2 A 

f 2 / 2 



U + + A j 
1 



J 



2 A-CZ + C/' + Aj + J 



(A14) 



The basis set is ^h+i> , |^et» ^o) , ^> , 

1^0, ^et) j l^et, ^L) , |^L, ^et)}- 

^Present address: Department of Physics, Hokkaido Uni- 
versity, Sapporo 060-0810, Japan. 
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